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ABSTRACT 

After defining geomagnetic noise as it applies to MAD, 
the geomagnetic indices currently used by the flee* to 
predict MAD geomagnetic noise are reviewed to determine 
their actual applicability. The current indices are deter- 
mined tc be insufficient, methods are proposed for estab- 
lishing a new MAD index, and a developmental MAD index 
system was tested. Geomagnetic fluctuations in the .04 to 
2.0 Hz frequency band were recorded at Monterey, California, 
and used for a preliminary test of the proposed MAD index. 
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I. GEOM AGNEr ICS REVIEW 
A. HISTORY OF GEOMAGNETICS 

The beginnings of the study of geomagnetism lie back in 
prehistory when magnetic attraction between iron and certain 
minerals was first observed. Exactly when this phenomenon 
was first noticed is not known, but the properties of magne- 
tite, then called lodestone, appeared in Greek literature 
around 600 3. C. (3rennan and Davis). [Ref. 1] 

Chapman [Ref. 2] indicates that the directional property 
of magnets was known and used in Europe prior to 1200 A. D. 
and possibly in China before then. E. N. Parker [Ref. 3] 
noted "It is an interesting fact that the ancient walls of 
Peking were lined up with magnetic north rather J -han 
geographic north, a difference at that time of about lo. We 
may presume That the surveyor found it easier to work with 
his compass needle by day than to sight on the pole star by 
nighT." This property allows the use of the Earth's geomag- 
netic field for navigational purposes. 

By the mid-fifteenth century it was determined in Europe 
that the magnetic compass does not point to true north. The 
angle between true north and the direction indicated by the 
compass is now known as magnetic declination by the geophy- 
sicist and as variation by the navigator. 

The magnetic field dip, or magnetic inclinarion is the 
angle, in a vertical plane, between the horizontal and the 
direction of the earth's magnetic field vector. It was 
observed in 1544 by an instrument maker in Nuremberg named 
Hartmann, and again by Robert Norman in London in 1561. 

These discoveries or observations gave rise to the study 
of geomagnetics as a specialty which gained its cornerstone 
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with William Gilbert in 160 0. After comparing his experi- 
mental results with the previous work of others such as 
Norman, Gilbert, in his book, "De Magnet e," concluded that 
"maqn us ma qr.es ipse est globus terrestris (the earth globe 
itself is a great magnet)’ 1 [Eef. 2]. It is this concept, 
that the earth is itself a magnet, that is the basis of the 
science of geomagnetics. 

Gilbert felt that the Earth's magnetism must remain 
constant except for geological changes, but it was soon 
determined that this was not the case. A ’secular varia- 
tion 1 of the Earth's was found to exist. 

Shorter term changes in the geomagnetic field were 
observed and it was eventually realized that geomagnetism is 
dynamic. In 1722 George Graham discovered that daily, or 
diurnal, variations exist [Ref. 2]. During the early 
Nineteenth century magnetic observatories began to be estab- 
lished to record the changes in the geomagnetic field in a 
systematic fashion (Knecht) [Ref. 4]. 

3. EARTH’S MAGNETIC FIELD 

1 • Co nstituent s of the Geomagnetic F ield 

There are various ways of breaking down the consti- 
tuting parts of the geomagnetic field. One way is no divide 
■he field in terms cf distance from the center of the earth. 
Doing this yields these three parts: internal, crustal, and 

external ( AFGL) [Ref. 5]. The internal field originates in 
the core region and is the more stable field, cottaininq 
only extremely low freguency temporal variations. The 
crustal (or anomalous) field arises from modifications made 
on the internal field by materials and structures in the 
Earth's crust. These variations are net spatially constant 
and give rise to some of what is known as geological 
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variations. The external field is the most dynamic and 
arises from many sources including the interaction between 
the solar wind and the terrestrial magnetic field. 

Another way of describing the components of the 



geomagnetic field is by time vari 


at ion . 


This divis 


ion is 


accomplished by considering that 


par 


t of 


the field 


wh ich 


varies with periodicities greater 


than 


about one year 


as t he 


steady field and what is left 


d 3 


t he 


variation 


field 


(Knecht) . [Ref. 4] 










The steady field con sis 


ts 


o f 


the above 


named 


internal field, also referred to 


d5 


the 


main field. 


Slow 


variations in the main field with 


peri 


ods 


of years or 


longer 



are referred to as secular variat ions . 
2. Models of the Sain Field 



Various models of the main geomagnetic field make 
use of a geocentric dipole. Gauss, in 1939, demonstrated 
that, as a fairly good first approximation to the geomag- 
netic field, the field of a uniformly magnetized sphere (for 
points outside the sphere) is equivalent to the field of a 
magnetic dipole located at the center of the sphere 
(Jacobs) . [ Ref. 6 ] 

The simplest of the present approximations of the 
geomagnetic field is that of a short bar magnet or dipole 
located at the center of the earth with an axis inclined 
approximately 11.5° from the Earth’s axis of rotation. The 
sense of the field lines is from south to north (Figure 
1 . 1 ) . 

The axis of this dipole intersects the earth at the 
geomagnetic north pole, 79. 5°N, 291. 0°E in geographic coor- 
dinates, and at the geomagnetic south pole, 78.5°S, 110. 0°E. 

The moment of the geomagnetic dipole is 9.1 x 10 22 amp-m 2 . 
It is these poles that are used to define the geomagnetic 
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Figure 1.1 Dipole appearance of Geomagnetic Field 

coordinate system (Knecht) [Ref. 4], The geomagnetic coor- 
dinate system is a spherical polar sysrem similar to the 
geographic coordinate system with a geomagnetic equator 
defined 90 degrees away from either geomagnetic pole in 
latitude. This tilted geocentric dipole model describes the 
main geomagnetic field to an accuracy of about 10%. 

In 1940 Chapman and 3artels defined an off-center 
dipole in the earth's interior, called the eccentric dipole. 
This dipole is displaced 0. 0635 earth radii (436 km) in 
magnitude from the center and in the direction of th=> point 
15. 6°N, 150. 9°S (geographic coordinates) (Vestine) [Ref. 7], 
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axis at the 



The intersections of the eccentric dipole 
earrh's surface are 81.D°N f 84. ana 75.3°S, 120. 40 e 

(Figure 1.2) (Haynes) [Bef. 8], This approximation is accu- 
rate to within a few percent. 




Pigure 1.2 Eccentric Dipole Model of Geomagnetic Field 



The field has additionally been modeled to an accu- 
racy of about 1^ by aeterai nining Gaussian coefficients by a 
least-squares fit of experimental measurements of “he 
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geomagnetic field. These coefficients are used in a spher- 
ical harmonic series representing the scalar potential of 
the field. This accuracy implies that the internal contri- 
bution to the total main field is at least on the order of 
99%. 

•The International Geomagnetic Reference Field (IGRF) 
yields values which differ by only parts per thousand from 
measured values. 

3. Sources of the Geom ague tic Field 

There are various elements that contribute to the 
geomagnetic field, some external to the earth's surface and 
some internal. As previously mentionned, the external 
contributions make up only a small fraction of the steady 
field, playing a more important role in the variation f ie ld. 

These external sources include current systems in 
the earrh's upper atmosphere affected by solar electromag- 
netic radiation and gravitation, solar corpuscular radiation 
or the interaction of solar plasma with the main field, and 
the effect of the solar interplanetary field. * Ref . 4] 

Various magnetic surveys of the world, including 
those conducted at ground level, by airborne instruments, 
and by satellite, have pointed to the fact that the largest 
source of the earth’s magnetic field is internal to it. 
While there exists residual permanent magnetism in the 
earth's crust, this cannot be the principal internal source 
of the geomagnetic field due to temperature and material 
properties known to exist in the earth's interior (Magata 
and Ozima) . [Ref. 9] 

Permanent magnetism is generated by microscopic 
electric currents, since a changing electric field will 
generate a magnetic field. Another way •‘■o generate a 
magnetic field is by the motion of electric charges in a 
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macroscopic current. Convective motion cf the electrically 
conducting fluid core of the earth, resulting in a macro- 
scopic current system, is considered to be the principal 
source of the main field. 

The most promising present theory of the generation 
of the geomagnetic field is that of some sort of a self-ex- 
citing dynamo system. This means that the motion of a 
conductor, such as the molten iron in the earth's core, in a 
magnetic field produces a current which in turn induces a 
magnetic field in support of the original magnetic field 
[Ref. 4, 7]. A very simple model of such a dynamo is shown 
in Figure 1.3 . 



O 




Figure 1.3 Simple Dish Dynamo 



The original excitation or seed field may be due tc 
an external field line, perhaps from the solar interplane- 
tary magnetic field. This original, poloidal, field line is 
wound up due to the differential rotation (notation not 
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constant with latitude) in the molten core. The wound up 
line becomes an intense azimuthal field which is carried 
outward by the upwelling associated with convection and 
twisted by the action of the Coriolis force. The twisting 
generates a helical toroidal field which, by outward diffu- 
sion, generates the externally observed quasidipole geomag- 
netic field. 

The combination and interaction of two (see Figure 
1.4) or more disk dynamos can also can also explain the 
reversal of the geomagnetic field [Ref. 7, 9]. 




con form 
iat ions 



Regional anomalies that are nondipolar (do not 
to the dipole field) possibly arise from eddy circu- 
in the outer core [ Hef . 4]. 
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• Magnetosp here 

The magnet osphe re can be defined to be that region 
(see Figure 1.5) occupied by the geomagnetic field above the 
ionosphere, a region where the field stongly influences the 
dynamics of ionized gas and charged particles (Kern) 

(Ref. 10]. 

If the space surrounding the earth were a perfect 
vacuum, the earth’s magnetosphere and magnetic field might 
be more or less symmetric and extend outward until it merged 
with, and its strength became insignificant compared to the 
solar and other planetary magnetic fields. This turns out 
not to be the case. 

Instead of being an island in a perfect vacuum, the 
earth encounters a continuous flow of hot, highly conductive 
ionized gas, or plasma, streaming outward from the sun. 
This continous stream of charged particles is called the 
sol ar wind. The density of this 'wind' near earth is on the 
order of 10 particles per cm 3 , and has a velocity averaging 
300-500 km per second (Jacobs). [Ref. 6] 

Both the solar wind and the geomagnetic field exert 
pressure. The hot plasma of the solar wind pushes against 
the geomagnetic field deforming the field. At distances 
greater than about 13 or 14 earth radii, the pressure of the 
solar wind greatly exceeds that of the geomagnetic field and 
the geomagnetic field will be swept along with the wind. 
From 8 to 10 earth radii inward, the geomagnetic field 

pressure will predominate, excluding the solar wind, this 
being the region of the mag netosp here. 

In the intermediate region, the magnitude of the 
solar wind and geomagnetic field pressures are comparable 
and the solar wind is compressed and flows around the 
geomagnetic field. This occurs when the magnetic energy 
density ahead of the plasma equals the kinetic energy 
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density of the streaming plasma. The solar wind is stopped 
at this point and forced to flow around the magnetosphere. 
This region where the magnetosphere starts is called the 
magnetop au se , [Eef. 8] 

The velocity of the undisturbed solar wind is analo- 
gous to a 'supersonic* velocity. Thus a shock front is 
formed between the magnetopause and the solar wind. The 
mag ne to s he ath is the region of severe turbulence that exists 
between the shock front and the magnetopause. 

Since the solar wind always travels outward from the 
sun, the affect of the wind on the earth's main field is not 
completely symmetric, although it is almost symmetric about 
an axis through the earth and sun. A ge oma gnetic tail is 
formed where the wind sweeps the geomagnetic field along 
with it cn the nightside of the earth [Ref. 4], Figure 1.5 
pictorially represents the effects of the solar wind on the 
geomagnetic field. 

5 . Ti me Va ri ations of the Geomagnat z.c Fie^d 

The geomagnetic field changes with time. As previ- 
ously mentioned, very slow variations in the main field with 
perioas cn the order of years to thousands of years are 
referred to as sec ular va r iati ons. Secular variations are 
geologic or ’ paleomagnet ic* in origin. Secular variations 
are not caused by a strength or orientation change of the 
center dipole. Paleomagne tic studies are used to determine 
the secular variation. Geologic structure, especially 
conductivity structure, may partially mask the secular vari- 
ation at one point on the earth as compared to that at 
another point. 

Other time variations of the field can be categor- 
ized into gu ie t variation fields and disturbed vari at ion 
fie Id s. Disturbed variation fields include geomagnet ic 
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Figure 1.5 Configuration of The Magnetosphere 

micro pal sat ions which will be discussed separately, since 
they are of particular in rerest as a noise source for MAD 
sensors. 



a. Quiet Variation Fields 

Quiet variation fields are those which are not 
due to disturbances in the interplanetary environment and 
which vary slowly and regularly [Hef. 4], 

There are several contributing fields to the 
quiet variation. These include the Solar Quiet Daily 
Variat i on (Sg) , the Lunar D ally Variation (L) , and the daily 
variation due to magn etosph eric effects. 
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The Solar Quiet (Sq) variation is the name given 
to the pattern of diurnal field variation with respect to 
solar local time which is caused by currents flowing in the 
ionosphere (Matsushita) [Ref. 11]. The major portion (about 
two-thirds) of the Sg field is due to what is referred as an 
atmospheri c dynamo. High speed tidal winds are generated by 
solar heating causing convection of the upper atmosphere 
[Ref. 4], These winds produce a stationary current system 
by moving the conducting particles of the upper atmosphere 
across the geomagnetic field lines. The daily variation is 
caused by the earth rotating under the current system. The 
remaining third of the Sq variation is caused by currents in 
the earth induced by the primary currents in the ionosphere. 

The Sq field can be shown to be latitude depen- 
dent reaching a maximum at the magnetic equator where a 
concentration of current, the equatorial electrojet, exists 
[Ref. 4]. The maximum horizontal component intensity is 
about 100 nT at the equator with 25 f o 50 n? more likely at 
higher latitudes. 

Longitudinal, seasonal, and solar cycle depen- 
dencies also occur for the Sq field. 

An example of the quiet-day variation at 
Monterey, California is shown in Figure 1.6 and is summa- 
rized in Table I. This data was taken using a Cesium Vapor 
total field magnetometer in February, 1979. 

The Lunar Daily Variation, L, is approximately 
one-tenth the magnitude of the Sq field and exhibits a 
semi-diurnal behavior in lunar time [Ref. 4], The major 
difference is that the winds are caused by lunar-solar 
gravitational tides. The L field is dependent on seasonal 
influences, lunar phase, the solar cycle, and latitude. 
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Figure 1.6 Variation of Total Geomagnetic Field Intensity, 
February 26, 1979 
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r ABLE I 



Quiet Day Variation in La Mesa Village, February, 26, 1979 

POSITION 36°36'N, 121°51’W (LA MESA VILLAGE) 

CALCULATED VALUES USING THE 1975 U.S. CHART MODEL (WORLO DATA 
CENTER A, BOULDER, COLORADO) FOR FEBRUARY 1979: 

0 I H(NT) ZCNT) F(NT) 



VALUE (FEB. 1979) 


: 15.96° 


60. 


74° 24, 


650 44,007 


50,44 


YEARLY CHANGE 


: -2.3 ' 


• • ; 


8 * -4 


t.3 -32.1 


-30. 


MEASURED VALUES 


OF TOTAL 


FIELD 


INTENSITY 


FOR FEB. 25, 


1979 


DATE 




TIME 




F(NT) 




i/25/79 




00:44 




50,405.24 








91:43 




50,405.79 








02:52 




50,404.80 








03:56 




50,405.48 








05:00 




50,408. 19 








06: OH 




50,409.15 








07: OS 




50,409.66 








01:22 




50,399.69 








09:26 




50,381.59 








10:20 




50,363.68 








2 1 : 2H 




50,345. 08 








12:21 




50,338.51 


MID-OAY l 






13:32 




50,356.48 








14:36 




50,378.86 








15:40 




50,386.98 








16:44 




50,384.30 








13:51 




50,400.81 








19:55 




50,396.91 








20:59 




50,401.35 








22:03 




50,400.30 








23:07 




50,400.54 




2/26/79 




00:21 




50,406.51 










AVERAGE 


50,390.50 ♦ 2 


1.43 (NT) 


f calc 


F M£AS = 


50.5 NT 


OR .11 
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A diurnal effect due to the dayside-n igktside 
difference in compression by the solar wind of the geomag- 
netic field causes a small variation of the order of 3 nT 

[Ref. 4]. 

b. Disturbed Variation Fields 

Disturbed varia tion fields are geomagnetic field 
variations that appear to be the result of interplanetary 
environmental changes and do not possess a simple period- 
icity. These variations include ionos pheric dis turban ces, 
the auro ra . geom ag netic st orms. and geomagnetic 
mic ro pulsations . 

An ion ospheric disturbance is a departure from 
the normal behavior of the ionosphere. 

Solar flare effects (SEE) are magnetic distur- 
bances produced by X-rays emitted from the solar flare. 
SFE’s usually have a rapid onset, typically a few minutes in 
duration, followed by a slower return to normal. The entire 
event lasts on the order of an hour (Reid). [Ref. 12] 

Auroras are caused by t he precipitation of 
charged particles down magnetic field lines into the atmos- 
phere and can be one of the brightest visual phenomena in 
the sky. The more intense and active auroras occur with 
geomagnetic disturbances and greatly increase icnizaticn as 
well as creating the spectactular visual displays. [Ref. 4] 

Geomagnetic storms are due to a change in the 
dynamic pressure of the solar wind. A typical storm begins 
with a compression of the magnetic field by an increase in 
solar wind dynamic pressure called a sudden commen ce ment 
(SC) , which increases the magnetic field (The so-called 
"gradual storm" begins with a gradual increase in field 
strength). The increase in field strength is on the order 
of several fens of nanoTesias (nT) and fakes about one to 
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six minutes to rise. If a disturbance starts with an SC but 
lacks the succeeding stages of a storm it is referred to as 
a s udde n im pu lse (SI) . Following the SC, the field remains 
compressed for two to eight hours in the initial ph ase of 
the storm. The main phase follows the initial phase. Over 
a period of hours to a day a westward ring current is set up 
at a distance of several earth radii whose magnetic field 
leads to a decrease in field strength on the order of 100 
nT. This decrease overshoots the eqoilibrium field strength 
and leads into the recover y phase of a day or longer where 
the field returns towards its prestorm strength as the ring 
currents gradually dissipate [Ref. 4], (Matsushita) 
[Ref. 13 ]. A magnitude-time graph of a typical geomagnetic 
storm is shown in Figure 1.7. 




Figure 1.7 Typical Sidlatitude Geomagnetic Storm 
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Geoaagn et ic »ic ropa lsations are rapid fluctua- 
tions in the surface magnetic field with periods of about 
0.2 seconds to 10 minutes (frequencies about 0.0016 to 5.0 
Hz) . These are observed as a type of geomagnetic distur- 
bance by ground based magnetometers [Bef. 4], 
Micropulsations will be discussed in depth later. 

6 • El ement s of t he Ma g neti c Field Vector 

The geomagnetic field vector is measured or charac- 
terized at any point by its direction and magnitude. This 
can be done in terms of some set of three independent param- 
eters such as two direction angles and the magnitude, or 
three perpendicular components. [Ref. 4] 

The system of coordinates commonly employed for 
describing the geomagnetic field on the surface of the earth 
is shown in Figure 1.8. The field is measured in terms of 
local (geodetic) coordinates with respect to True North. 

The various coordinates are referred to as magne tic 
ele ment s and are defined as fellows: 

B: Total Field Intensity (the symbol F is 

also used) 

H: Horizontal Component 

X: Northward, or North-South Component 

Y: Eastward, or East-West Component 

Z: Downward, or Vertical Component 

D: Declination or magnetic variation. 

This is the angle between X and H and is 
measured positive eastward. 

I: Inclination or dip angle. This is the 

angle between H and 8 and is measured 
positive downward. 
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Figure 1.8 



Haguetic Field Elements 
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II. INTRODUCTION TO HAG NE TIC ANOMALY DETEC TION (MAD) 

A. DEFINITION OF A MAGNETIC ANOMALY 

A magnetic anomaly is defined as any spatial variation 
cr disturbance in the geomagnetic field which is due to 
local causes. Anomalies can be caused by waves, ere 
deposits, sea mounts, and magnetized objects such as surface 
ships and submarines (Anderson) [Ref. 14], For the purposes 
of this research, magnetic anomalies due to a simple dipole 
field, such as those generated by submarines, will be 
regarded as signals, while other anomalies will be regarded 
as noise and will be discussed later. 

3. HISTORICAL DEVELOPMENT OF THE MAGNETIC DETECTION OF 

SUBMARINES 

1 • D etection Sys tems 

Attempts at finding submerged submarines by sensing 
disturbances in the geomagnetic field date bach as least as 
far as World War I. In 1918, Earnest Merritt, at the Naval 
Experimental Station, New London studied the use of a fixed 
coil type of detector for use in moving boats and airplanes 
(Slichter) . [Ref. 15] 

MAD, originally hr.own as Magnetic Airborne Detector, 
and now called Magnetic Anomaly Detection, had its begin- 
nings as an airborne ASW sensor in late 1940 and early 1941. 

Anomalies in the geomagnetic field caused by the 
presence of submarines, are on the order of one to a few nT 
(gammas) in magnitude. This is quite small compared to Mi? 
magnitude of the field itself (30000 - 60000 nT) . Two 

methods are generally employed to measure this small 
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disturbance. One is to use a gradiometer which measures the 
spatial rate of change of the magnetic field or its 
gra dient . The second method, and the one presently used by 
U.5. Navy aircraft is to use a magnet ome ter to directly 
measure changes in the magnitude of the magnetic field. 

The British investigated the use cf a gradicmeter 
fcr submarine detection, and by early 1941, had developed a 
two-ccil gradiometer system which could detect a submarine 
at a range of 200 fee-, under favorable conditions. This 
range was considered to be too small to be of operational 
value and work, on such a system was terminated when a magne- 
tometer system showed premise (Coleman). [Ref. 16] 

By late 1940 , Victor V. Vacquier of the Gulf 
Research and Development Company had developed a sensitive 
sat ur a bl e core magnetometer intended for geophysical 
(mineral) prospecting. The Vacquier magnetometer became the 
basis for further MAD development. The Airborne Instrument 
Laboratory of Columbia University continued the investiga- 
tion of means of localizing submerged submarines by MAD. 

The simplest saturable core or fluxgat? magnetometer 
consists of a saturable or ferromagnetic core around which a 
coil of wire is wrapped (see Figure 2.1). This coil carries 
a sinusoidal current, I (t) , which is large enough to satu- 
rate the core during part of each cycle. The inductance of 
the coil will change as a function of the magnetization of 
the core. The core magnetization, in turn, depends on the 
instantaneous current in the coil, and, if present, the 
external magnetic field. 

In the absence of an externally applied magnetic 
field, magnetization as a function of exciting current is 
symmetric around 1=0 (Figure 2.2). An external magnetic 
field parallel to the core's axis will change the magnetiza- 
tion of the core and shift the magnetization curve. This 
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Figure 2. 1 Element of Fluxgate Magnetometer 

change will cause an asymmetry which can be sensed by 
analyzing the harmonic content of the signal. The coil can 
be combined with a stabilizer system, which keeps the 
detector element (coil) aligned with the geomagnetic field. 

This was the basis of the MAD magnetometers used in 
World War II such as the AN/ASQ-1, ASQ- 1 A, and ASQ-2 

[Ref. 16]. 

The fluxgate magnetometer measures only the compo- 
nent of the external field parallel to the axis of the 
ferromagnetic core. In order to measure the total field in 
this fashion if is necessary to align the ferromagnetic core 
along the earth's magnetic field or by using mutually 
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Figure 2.2 Magnetization Versus Exciting Current 

perpendicular (orthogonal) cores. Precision requirements 
made this type of system difficult to realize during World 
War II but it has since found many years of operational 
usage and has found widespread use in geophysical explora- 
tion work and satellite mapping of the geomagnetic field. 

2 • 5a r ly Operation al 3 sage 

During World War II MAD provided a passive raerhcd of 
detection and tracking of submerged submarines. Within its 
range (then on the order of 500 feet) MAD gave a measure of 
surprise to the attacking aircrew. Until the first attack 
was delivered the crew of the submerged submarine might no* 
even be aware of the aircraft hunting it. 

Operational deployment of MAD began on a limited 
basis in December, 1941 with the installation of the early 
Mark I MAD in a blimp at Naval Air Station Lakehursr, New 
Jersey [ Ref. 16 ]. 
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3y the end of 1942, MAD was operational on beard PBY 
Catalina aircraft, nicknamed "Madcats. " Even though MAD was 
operational at this time, it was not until February, 1944 
that an initial contact by MAD led to the sinking of a 
sub marine. 

During February, 1944, the Hadcats of Patrol 
Squadron 63 were assigned to fly a MAD barrier patrol of the 
Straits of Gibraltar. On 24 February 1944, a Madcat of 
VP-63 detected U-761 by use of MAD and commenced tracking to 
confirm the contact as a moving target. An attack was 
conducted in conjunction with another Catalina, two 
destroyers, and eventually two other aircraft, and U-761 was 
sunk. U-392 was sunk after a similar MAD contact on 16 
March, and on 5 May, the third successful attack on a U-boat 
resulting from an initial MAD contact took place when U-731 
was sunk (OEG No. 51, Price). (Ref. 17, 18] 

The attacks at the Straits of Gibraltar demonstrated 
that though MAD has a limited search rate, there are 
scenarios where it can be employed effectively as a search 
sensor, such as providing a blockade across a restricted 
area without the presence of surface craft (DEG No. 54) 
(Ref. 19]. 

3 ♦ Current S vs terns 

The magnetometer system in current operational use 
is the optically pumped magnetometer. The optically pumped 
magnetometer measures the external magnetic field by making 
use of the fact that when an atom is immersed in a magnetic 
field, its energy levels are split. This is known as the 
Zeeman effect. For the fields of interest, the amount of 
splitting of the levels is proportional to the intensity of 
the magnetic field. By measuring the separation be i !w°en the 
levels the magnitude of the magnetic field can be 
det er mined . 
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This type of magnetometer usually makes use of 
Cesium or Rubidium vapors, or Helium gas. The current oper- 
ational MAD system, the AN/ASQ-81, is a Helium gas optically 
pumped magnetometer. 

In an optically pumped magnetometer (see Figur° 
2.3), the sample vapor or gas (such as Helium) is collected 
into an absorption cell. Circularly polarized light is 
passed through the cell giving up some of ins energy in 
exciting or pumping the electrons of the sample gas to 
higher energy levels. These electrons then fluoresce to 
lower, metastable energy states. This is 'optical pumping.' 
A detector monitors the degree of optical pumping by meas- 
uring the transparency of the gas cell. 




Figure 2.3 Metastable Helium Magnetometer 



The actual separation of the energy levels is meas- 
ured by applying a weak R. F. field which redistributes the 
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electrons among the ground stare sublevels. These tran- 
sitions will only occur when the R. F. field has a partic- 
ular frequency (called the Larmor frequency) which is 
proportional to the separation between the levels and there- 
fore to the external magnetic field [Ref. 20]. Thus meas- 
urement of this frequency yields a direct value of the 
magnetic field. 

4 . Fut u re Developments 

Optically pumped magnetometers are the most sensi- 
tive sensor for fleer usage today. Increasing capability to 
make use of the phenomenon of superconductiviry will yield 
more sensitive MAD sensors. Superconducting Quantum 
Interference Devices (SQUIDS) which make use of the 
Josephson effect have a theoretical sensitivity of 10~ 7 nT 
as opposed to the theoretical value of 0.001 r.T and the 
operationally realized value of 0.01 nT in current fleer 
systems (Chilton) [Ref. 20]. Field changes on the order of 
10~ s nT have actually been measured. This SQUIDS are used 
in both superconducting magnetometers, and superconducting 
gra diomete rs. 

A Josephson junction consists of a thin layer of 
insulator between two superconductors. This junction has 
the property that a current can flow across it without 
developing a voltage up to some maximum current. A voltage 
is developed for all current values greater than the maximum 
current. 

In a superconducting magnetometer a superconducting 
ring containing a pair of Josephson junctions is used to 
measure the amount of magnetic flux penetrating the loop. 
This system is a vector magnetometer which measures varia- 
tions in only one field component. 
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3y using more than one superconducting loop, it is 
possible to measure magnetic field gradients, and construct 
a superconduc ting gradiometer. This is also a component, or 
vector, sensor. 

These instruments have been used extensively in the 
laboratory and to some extent in geophysical work, but 
considerable engineering problems remain to be solved before 
operational Navy use can be contemplated. 

C. MAD SIGNAL AND BANDPASS 
1 . Source of the Signa 1 

The NAD signal results from moving a magnetometer 
through the magnetic field of a submarine, which can be 
approximated by the field of a magnetic dipole. Figure 2.4 
depicts the formation of a submarine caused anomaly. 

The magnetic moment of an object in the earth's 
magnetic field can be due to permanent magnetization, 
magnetization induced by the earth's magnetic field, or a 
combination of bcth. In the case of a submarine hull both 
causes are present with a small amount of permanent 
magnetization produced by hull stress in metal components 
during construction and stress caused bv submarine diving 
and surfacing. This, however, is a minor contribution. The 
most important constituent of the submarine magnetic moment 
is due to magnetization induced in the hull by the presence 
of the geomagnetic field. [Ref. 21] 

The induced field of a submarine depends on the 
effective permeability along the vertical, athwartships , and 
longitudinal axes of the submarine. This information , 
taken together with the strength and dip ar.ale of th=> 
geomagnetic field, yields a fairly precise calculation of 
the magnetic anomaly produced by ohe submarine. ’ Depermir.g' 



36 




Figure 2.4 Formation of Submarine Caused Magnetic Anomaly 
(Simplified 7iew Sear Equator) 
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of the submarine hull cancels out the permanent 
magnetization leaving the induced magnetization as the 
principal signal source. [Ref. 21] 

The MAD signal is approximated as the projection of 
the submarine dipole field onto the geomagnetic field 
vector. This approximation is good because the magnitude of 
the dipole field is very much smaller than the magnitude of 
the geomagnetic field (a few nanoTeslas as opposed to 
approximately 50000 nT) . Therefore, whenever the dipole 
field is perpendicular to the earth's field a region of zero 
signal will result. The signal recorded by the AN/ASQ-81 or 
other MAD equipment is a mapping along the aircraft's flight 
path of this signal. Figure 2.5 qualitatively describes 
some aspects of the MAD signal. 

2 • An d er so n Functi ons 

The submarine anomaly signal shape is a function of 
the dip angle of the geomagnetic field, the magnetic heading 
of the aircraft, magnetic heading of the submarine dipole, 
and the lateral range between the aircraft and the subma- 
rine. These factors detwermine the 'A' coefficients for the 
Anderson functions below. 

In 1949, J. E. Anderson of the Naval Air Development 
Center determined that the MAD signal, obtained along any 
course, consisted cf a linear combination of three basic 
components. Different shaped signals could by changing the 
proportions! contribution of these basic components. The 
mathematical representation of these components are now 
referred to as the Anderson functions. Using the dimension- 
less parameter 'b* (defined as the distance traveled along 
the aircraft trach divided by the slant range at closest 
point of approach (CPA) , b=0 at CPA) the anomaly can be 

represented by 
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Figure 2.5 Qualitative Aspects of MAD Signal 
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M A2 b 2 + A 1 b + AO 

Bs (2- la) 

5/2 

(b 2 +1) 

or 3S [AO £0 (b) * A, f1(b, * A2 «(M] <2-1b» 

where M = magnetic moment of target diode 

Z = lateral range between target* and plane at CPA 

f 0 (b) = 

5/2 

(b 2 ♦ 1) 
b 

f 1 (b) = 

5/2 

(b 2 ♦ 1) 
fc 2 

f 2 (b) = 

5/2 

(b 2 ♦ 1) 



The strength of the detected signal is seen to fall 
off as the cube of the distance between the aircraft and the 
target. 

By analyzing the functions and their coefficients, 
it can be shown that the optimum orientation for maximum 
anomaly signal defection is when the aircraft, tarqet 
dipcle, and geomagnetic field are lined up together as 
closely as possible. Specifically, this occurs when the 
submarine moment and the aircraft's track are oriented 
North-South. [Ref. 14, 21] 

The Fourier transform of the component functions 
were taken to determine the frequency distribution of energy 
in the a AD signal. The anomaly signal components are shewn 
in Figures 2.6 through 2.8 (Anderson) [Ref. 14]. The 
Fourier transforms of the signal components with a platform 
velocity of 150 knots are shown in Figure 2.9. 
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Figure 2-6 Basic MAD Component fO (b) 
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Figure 2-7 Basic HAD Component f1{b) 
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Figure 2.8 Basic MAD Component f 2 (b) 
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Figure 2.9 Frequency Spectrum (Fourier Transform) of MAD 
Signal Components 
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The MAD Filter Ban d pass 

The MAD signal, being a transient, is essentially a 
broadband (as opposed to a discrete frequency) signal. The 
frequency spectra of energy content shows how the signal is 
contained in a certain frequency range. In order to screen 
out unwanted noise, bandpass filtering is used in processing 
magnetic anomaly detections. 

The major factors influencing the determination of 
an optimum MAD filter frequency are aircraft speed and the 
slant range from the aircraft to the target at CPA. Factors 
exerting a minor influence include dip angle, aircraft 
heading, and target dipole orientation. 

The range in frequency variation is about an octave, 
with the higher values occuring when the passes are made 
parallel to the axis of the dipole and the lower when the 
pass is perpendicular to the dipole. Anderson empirically 
determined that the center frequency of the filter is given 
by 

0.4 v 

f = (2-2) 

Z 

where f = frequency in Hertz 

v = aircraft" velocity in ft/sec 
and Z = greatest anticipated CPA 
range in feet 

The ASQ-81 3andpass Filter has highpass settings of 
0.04, 0.06, 0.08, and 0.1 Hz, and lowpass settings of 0.2, 
0.4, 0.6, and 2.0 Hz (Orion Service Digest 26) [Ref. 22]. 

There are no recommended settings for normal operation since 
background noise varies. The filter characteristics for the 
0.06 to 0.6 Hz settings are shown in Figure 2.10 (Orion 
Service Digest 28) [Ref. 23]. 
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The figure is a representation of the adjustable band pass 
filter of the ASQ-81 by itself. Other parts of the ASQ-81 
system add in an additional high pass filter which adds 
another 12 dB/octave roll off to the low frequency end in 
Figure 2.10. Thus Figure 2.10 would represent the charac- 
teristics of the entire ASQ-81 system if the roll off of 36 
dB/octave is changed to 48 dB/octave. [Ref. 24] 




Figure 2.10 ASQ-81 Filter Characteristics: 0.06 - 0.6 Hz 

Using equation 2-2, the center frequency for a MAD 
signal with a CPA slant range of 300 feet and aircraft 
velocity of 220 knots is 0.49 Hz. Using a lcwpass filter 
set ring above 0.6 Hz would not be very useful. 

One reason the ASQ-81 filter extends up ro 2.0 Hz 
instead of just 0.6 Hz is the following. The early Service 
Test Engineering Model (STEM) ASQ-81 that preceeded the 
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current ASQ-81 production model had two bandpass filters, 
one adjustable from .04 to .6 Hz, and the ether fixed at .75 
to 10 Hz. During tests, the .75 to 10 Hz channel proved to 
be helpful in monitoring the STEM ASQ-81 "system" noise 
level. During design reviews for the production ASQ-81, it 
was decided that it might be useful to retain the ability to 
monitor "system" noise. In order to do this easily, -he .75 
to 10 Hz band was dropped and a 2 Hz position was added to 
the adjustable filter. [Bef. 24] 

The center frequency at 1200 feet and 160 knots is 
.09 Hz. Consequently, the highpass filter settings could 
affect slower, longer range MAD detections. 

The choice of the highpass (lower end of bandpass) 
settings have to be made after considering the noise present 
at time of system operation. 
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III. SOURCES OF MAD NOISE 



A. INTRODUCTION 



Or. a' practical level MAD noise is defined as magnetic 



disturbances falling within the MAD passband (0.04 to 2.0 
Hz f or since 0.6 Hz is the normal upper limit, 0.04 to 0.6 
Hz) and having an amplitude greater than 0.01 nanoTeslas 
(the system sensitivity of the AN/ASQ-31 MAD system) . 

MAD noise sources can be divided info the following 
categories : 

-Equipment Noise 
-Aircraft Platform Noise 
-Aircraft Manuever Noise 

-Gradient Noise Due to Aircraft Motion Through the 
Geomagnetic Field 
-Geologic Noise 

-Noise from Hind Waves and Swells 
-Geomagnetic Noise 

3. SENSOR, PLATFORM AND MANUEVER NOISE 



Sensor noise is the self-noise generated by the opera- 
tion of the equipment itself. This can be partially due to 
the fact that the detect or element is misaligned with 
respect to the geomagnetic field vector. Changes in lamp 
intensify, phofodetector noise and noise in the electronic 
circuits can also contribute to sensor noise [Ref. 21]. The 
self-noise limitation of the ASQ-81 is 0.01 nT. 

Platform noise is generated by components fixed to the 
aircraft in the vicinity of the sensor [Ref. 23]. 
Permanent, induced and eddy-current magnetic fields are 
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associated with the airframe. Permanent magnetic fields are 
due to aircraft structure or equipment having ferromagnetic 
parts. This field changes its orientation with respect to 
the geomagnetic field vector as the aircraft manuevers, 
causing field fluctuations near the magnetometer. 

P la t form noise is also induced in aircraft ferromagnetic 
structures by the geomagnetic field. Similarly, eddy-cur- 
rents are induced in aircraft skin, ribs, and frames, and 
these currents , in turn, cause additional magnetic fields. 
Thus, rapid aircraft manuevers will induce changes in the 
magnetic field sensed by the magnetometer. 

Platform noise in aircraft mounted sensors is countered 
by applying equal and opposite magnetic fields to the sensor 
in a process called compensation. Towed MAD systems are 
essentially free from this type of noise. 

C. GRADIENT NOISE 

Gradient noise car. be divided into turn noise and noise 
due to changes in altitude. 

Turn noise is a problem when 'MAD trapping' or using MAD 
for tracking a target. The earth's magnetic field has a 
horizontal gradient (in this case the magnitude varying with 
latitude). As the aircraft moves in the direction of the 
gradient the field strength changes. The noise due to a MAD 
trapping or hunting circle is centered in frequency at the 
reciprocal of the time taken to complete or.e revolution. In 
the case of a two-minute circle, the noise would be centered 
at 0.00833 Hz, well below the filter used in the ASQ-81. 
The horizontal gradient noise due to flying a clover leaf 
pattern, for the most part, also falls below the MAD 
passband. [Ref. 21] 

Vertical gradient, noise is due to changes in sensor 
altitude. An altitude gr a die no of up to 0.005 to 0.01 nT 
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In arsas of 



per foot exists in the earth's main field, 
geological anomalies this gradient is even larger. Fast 
altitude or aircraft pitch changes can cause a magnetic 
field fluctuation of sufficient amplitude to be of concern. 
To avoid vertical gradient noise, altitude compensation 
equipment is used. [Bef. 23] 



D. ENVIRONMENTAL NOISE 

Magnetic noise from sources existing in the natural 
environment include geologic noise, temporal variaticn in 
the earth's magnetic field, and noise due to ocean waves and 
swe 11s. 

Geologic noise has its source in naturally occuring 
magnetic anomalies caused by magnetic material present in 
the earth's crust. When th® sensor passes over geological 
anomalies, the relative motion causes a MAD-lik® signal to 
be recorded. Geologic noise is usually more pronounced in 
shallow water as the sensor is much closer to the source of 
the noise. Geological magnetic anomalies are often associ- 
ated with such oceanographic features as seamounts and ocean 
rid ges. 

Sea water is a conducting medium which is transported by 
the physical motion of water waves in ^he presence of the 
geomagnetic field. This motion induces currents in the sea. 
These currents give rise to secondary magnetic fields, which 
add vectorially to the quasistatic, geomagnetic field 
(Weaver) [Ref. 25]. These fields can be detected at signif- 
icant distances above the sea surface and fall off exponen- 
tially with altitude. Figure 3. 1 is a plot for several 
surface wave periods of the induced magnetic field per meter 
amplitude of the surface wave. These induced fields can be 
a problem at the lew altitudes where MAD is used. 
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Figure 3.1 Induced Magnetic Field per Meter of Saveheight 
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E. GEOMAGNETIC NOISE 

Temporal variations in the earth’s magnetic field with 
frequencies in the MAD bandpass and amplitudes greater than 
0.0 1 nT have become known as g e omagnetic noise in the MAD 
literature. 

Quiet daily variations, such as the Sq and L variations, 
have periods sufficiently long to fall far below the MAD 
passband . 

Geomagnetic storms have been discussed previously. 
Sapid fluctuations with high amplitude and falling within 
the MAD passband occur in connection with geomagnetic 
storms. 

G eomagr et ic m icropulsa t ions comprise the last category 
of geomagnetic noise to be discussed. 

1 . Ge omagnet ic Mien ogu Is at ions 

Geomagnetic micropulsations are rapid fluctuations 
of the earth's magnetic field with periods from 0.2 seconds 
to 10 minutes and amplitudes from about 0.1 nT to as high as 
a few tens of nT's. These fluctuations are caused by elec- 
tromagnetic perturbations propagating in the magnetosphere 
as hydromagnetic waves (Nishida). [Ref. 26] 

Micropulsat ions are classified by morphology, that 
is, by examining periods, amplitudes, times of occurence and 
other observed charac terist ics. Micropulsations have been 
placed into two broad categories: ir regular and c on ti nuo us. 

The irregular pulsations are represented by the symbol Pi, 
and the continuous pulsations by Pc. Table II shows the 
breakdown by period of the various micropulsation classes. 

It can be seen that for MAD geomagnetic noise the 
lower frequency Pci, Pc2, Pc3, and Pil pulsations are of 
interest . 
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TABLE II 



Geomagnetic Nicropulsation Classes 



NOTATION 


PERIOD (sec) 


AVERAGE AflPLITUD 


Pci 


0.2 - 5.0 


0.05 - 0.1 


Pc2 


5 - 10 


0.1 - 1.0 


Pc 3 


10 - 45 


0.1 - 1.0 


Pc 4 


45 - 150 


0.1 - 1.0 


Pc 5 


150 - 60 0 


1 - 10 


Pil 


1 - 40 


0.05 - 0.1 


Pi2 


40 - 150 


1 - 5 



(after Jacobs) 

Pci pulsations are regular sinusoidal oscillations 
with periods nornallly falling in the 0.3 to 4 second (0.25 
to 3.33 Hz frequency) range. Pci’s may start as separate 
bursts and gradually develop into a series of pulsations 
which could last for hours. They may also occur as consecu- 
tive groups of pulsations with sharply varying frequencies. 
The average amplitude of Pci’s is 0.05 to 0.1 nT and they 
tend to have a single well defined frequency. Pci pulsa- 
tions with frequency less than 0.5 Hz are more common at 
high latitudes than at mid latitudes. These pulsations 
occur in the daytime in the auroral zone, and at night and 
early morning hours in lower latitudes. Pci’s are charac- 
teristic of the quiet and weakly disturbed states of the 
geomagnetic field and show an increase in activity one to 
two hours before, and four to seven days after, a magnetic 
storm. (Ref. 6] 

Pc2 and Pc3 pulsations are grouped together in most 
characteristics. The amplitudes of the oscillations are 
usually under 0.5 nT and the typical frequency range is 0.03 
to 0.2 Hz. These pulsations are normally found during the 
day with their activity reaching a maximum around noon. Pc3 
pulsations show a seasonal variation with a minimum of 
activity occuring during winter. 
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Pc4 and Pc5 are large amplitude f luctuarions , but 
fall below the frequency band of interest for MAD. 

Pil pulsations have an irregular form with an 
average amplitude of 0.01 to 0 . 1 nT and a frequency mainly 
in the 0.10 to 0.17 Hz range. Spectral analysis of these 
pulsations show a braod band of frequencies. Pil amplitudes 
have maximum values in the auroral zones, with the intensity 
of the pulsations decreasing with decreasing latitude. 
Pil's are normally observed in the late night and early 
morning hours, and show an increase in activity with 
increased geomagnetic field disturbance. [Ref. 6] 

Geomagnetic micropulsations can be observed anywhere 
on the globe, at various times of day and year, and in both 
quiet and disturbed geomagnetic field conditions. 
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IV. 



METHODS OF EVALUATING GEOMAGNETIC ACTIVITY 



A. 



I NT ROD OCT ION 



Previous chapters have defined geomagnetic activity as 
it applies to Magnetic Anomaly Detection. In this chapter, 
methods of evaluating that activity will be examined, 
including methods currently in use by the fleet. 



B. GEOMAGNETIC INDICES 

A geomagnetic index is simply a measure used to quantify 
and describe time variations of the earth's magnetic field 
resulting from solar-terrestrial relationships. These 
indices are commonly used to express the intensify and 
depict the character of geomagnetic activity throughout the 
day . 

For the most part, geomagnetic indices developed as 
range indices, measuring the difference be tween the high and 
low values for different field components measured during 
the day by magnetic observatories (Lincoln) [Ref. 27]. Most 
current indices are of the range type, but other indices 
have been developed which are more subjective or qualitative 
in nature. 

Geomagnetic indices are designated by a lerter coda such 
as: C, Ci, Cp, C9, Q , R, w measure, Dst, K , Ks, Kp, ak, Ak, 

ap, and Ap. There are add! tonal indices in use. 

The C, Ci, Cp, and C9 indices are daily magnetic field 
character figures. The C index is the daily character 
figure for a single observatory. In this scale, C=0 indi- 
cates a quiet day, C= 1 a moderately disturbed day, and C=2 a 
heavily disturbed day. The daily international character 
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figure, Ci , is the arithmetic mean of the C indices reported 
by participating observatories around the world. Cp, the 
daily planetary character figure, is similar to Ci, except 
that it is derived from the values of Kp and ap. C9 is a 
contracted scale for Ci an Cp with single digit values 
running from 0 to 9 (Bartels). [Ref. 28]. 

The Q and R indices are quarter-hourly and hourly range 
indices respectively, taken at high latitude stations only. 

The W measure is an index of the equatorial electrojet. 
Dst is a measure of ring current effect. They are both 
amplitude indices. [Ref. 27] 

1 . K, a, and A I ndi ces 

The K, Ks, Kp , ak, Ak, ap, and Ap indices comprise a 
group of related 3-hour range indices. The K index is a 
single station code using a quasi-logarithmic scale from 0 
to 9 to measure geomagnetic activity. The value of K is 
determined by first determining the difference between the 
lowest and highest deviations from the regular daily varia- 
tion (Sg) during a 3-hour period. This range (in nT) is 
converted to the K scale based on the historical activity 
ranges at the particular observatory involved [Ref. 27]. 
The conversion for the Fredericksburg, Virginia observatory 
is given in Table III. This conversion can also be applied 
to the USAF/NOAA observatory in 3culier, Colorado [Ref. 29], 

The Ks index is a standardized K index which is 
freed from local variations and is then used to determine 
the planetary 3-hour index, Kp. 

The equivalent three- hour-range, ak, is a conversion 
of the K index as shown in Table IV. In order to detemine 
the units of ak for a particular observatory, divide the 
lower range limit of K=9 by 259. Thus for Fredericksburg 
and Boulder, ak is in 2-nT units. 



TA3L2 III 



Conversion from Bangs to K for Fredericksburg, 7a. 



K 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 



Ra r.ge^nl) 



5 - 
10 - 
20 - 
40 - 
70 - 
120 - 
200 - 
330 - 
>500 



10 

20 

40 

70 

120 

200 

330 

500 



(after Lincoln) 



TABLE IV 



Equivalent Bangs ak for Given K 



K0 123 456 78 9 

ak 0 3 7 15 27 48 80 140 240 400 

(after Lincoln) 



Ak is the equivalent daily amplitude and is the 
average of the eight daily ak valuas at a particular observ- 
atory. This index is promulgated using the name of ths 
observatory, the Ak index for Fredericksburg is known as the 
A-? rede ricks burg or A-Fred index. 

The equivalent planetary amplitude, ap, is deter- 
mined from the Kp index in a fashion similar to chat of 
determining ak from the K indices. The eight ap values for 
a given day can then be averaged into the daily equivalent 
planetary amplitude A p. These two indices are given in 2-nT 
units . 
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C. GEOMAGNETIC INDICES IN FLEET USE 
1 • Cu rrent Usage 

Geomagnetic indices of interest, in connection with 
MAD operations are the K, A k, and Ap indices. 

Fleet operators utilize the " alpha Index' for 
perdicting geomagnetic activity over the entire world 
(Ref. 30]. This index is promulagated by the Fleet 
Numerical Ocenaograph ic Canter, Monterey, California, in the 
environmental briefings received by aircrew personnel. This 
index is the Ap index as sent out from the Space 
Environmental Services Center, Boulder, Colorado in the 
Joint USAF/NOAA Report of Solar and Geophysical Activity 
[Ref. 31]. 

The Boulder K index is available to interested 



parties 


by 


telephone 


recor 


ding and 


in 


the WWV 


and 


WW V H 


radio 


bro; 


ideas 


ts 


and is 


there 


fore aval 


lable t 


o f 


leet 


users 


[Re: 


f. 32 


]• 






















The 


Ak inde 


x fo 


r the 


Fra 


dericksburg , 


Vir 


ginia 


cbs ; 


er vat 


cry 


has been 


used 


in studi 


es 


of geema 


gnat 


ic act 


ivitv 


as < 


appli 


ed 


to Magneti 


c Ano 


maly Det 


act 


ion [ Ref 


. 33 


]• 






2. 


The 


oret ical Applic 


a bili ty 


2? 


A and K 


Indi 


ces to 


MAD 






With the wid 


est 


u se fu 1 f 


lit 


er sett 


ings 


, the 


MAD 


ban < 


ipass 


r a 


nges from 


0. 04 


to 0.6 


Hz 


(1.7 to 


25 


seconds in 


per: 


Lod) . 




As such. 


in or 


der for 


a 


geomagne 


tic 


index 


to be 



directly applicable for MAD use, it should be sensitive to 
that frequency range. 

The K indices, and the K-derived A indices, are not 
especially sensitive to the MAD range. Mayaud [Ref. 34] 
indicates that these indices are mainly sensitive to fluctu- 
ations whose periods are much longer than the lower end of 
the frequency range analyzed, that is, a frequency corre- 
sponding to a period of 45 minutes (0.0004 Hz). 
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One reason for this lack of sensitivity for MAD 
bandpass geomagnetic noise is that the amplitude of geomag- 
netic fluctuations varies inversely with frequency, so that 
the amplitude of the fluctuation increases as the as the 
frequency decreases. If can therefore be seen that the 
fluctuations with periods of an hour or greater largely 
determine the variation range used to calculate the K index. 
The activity driving the K and A indices, will, because of 
band pass filtering, not even be observed by the MAD system, 
and the activity of interest to MAD might not influence the 
K or A indices at all. It can be concluded that there is no 
direct physical link between MAD geomagnetic noise ana 
either the K or A indices. [Ref. 35] 

3. Exp erimen tal Correlation of A and K Ind ices with MAD 
Band Noise 

a. ASQ-10A Study 

Brennan and Smifs [Ref. 33] found that geomag- 
netic micropulsation activity was recorded at their ASQ-10A 
MAD magnetometer site in Maryland, when the A-Fredericksburg 
index was greater than 25. This cccured everytime they were 
recording data with A -Fred greater than 25. It is important 
to note that they observed additional activity during some 
periods when A- Fred was lass than 25. 

The 3renrar. and Smits study tends to validate 
the use of the A indices as at least qualitative indications 
of geomagnetic noise in the MAD bandpass. Their study, 
however, was specific to the ASQ-10A magnetometer system 
which has a sensitivity of 0.1 nT as opposed to the 0.01 nT 
sensitivity of the ASQ-81 system. The effect of the ASQ-10 
sensitivity is to filter out most Pci pulsations. As Pel 
pulsations do not correlate well with the A and K indices, 
the filtering out of these pulsations would tend to increase 
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the reliability of the A-Fred or Ap indices as measures of 
HAD geomagnetic noise. Mason [3ef. 36] stated that the 
"occurence of Pci is well known to be associated with low Kp 
values." It should be noted that an operational drawback of 
the filtering out of Pci pulsations by the ASQ-10A is than 
the system also filters out valid signals of less than 0.1 
nT amplitude. 

While information has been presented that 
suggests the the A-Fred index can be useful for MAD geomag- 
netic noise evaluation for the AN/AS2“10A system, sufficient 
data was not presented to draw conclusions for index usage 
with the AN/ASQ-8 1 system. 

b. ASQ-81 Study 

For two weeks in April 1976, Naval Air 
Development Center personnel operated a geomagnetic observa- 
tory at the Atlantic Undersea Test and Hvaluation Center in 
the Bahamas. The primary magnetometer used for this observ- 
atory was the ASQ-81 magnetometer. One purpose of this 
observatory was to compare the K index (as determined by the 
San Juan Observatory) with geomagnetic activity in the HAD 
band pass. The conclusion of this study was that the K-San 
Juan index did not correlate with geomagnetic noise in the 
MAD band pass. [Ref. 37] 

We have undertaken a correlation analysis using 
the Fisher Z transformation [Ref. 38] of the NADC data 
furnished to us by Ochadlick. The data used is listed in 
Tables V and VI. 
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TABLE V 



Observed 


AUTEC Data 


and Indices, 


April 


11-18, 1976. 




Time (Z) 
03 00 






Maximum 


Date 

1 1 


K-San Juan 
1 


ts 


Amplitude (nT) 

0.17 


1 1 


0600 


2 


12 


0.14 


1 1 


1500 


2 


12 


0.62 


12 


00 00 


2 


9 


0.21 


12 


03 00 


3 


9 


0.28 


12 


06 00 


3 


9 


0.28 


12 


1500 


0 


9 


0.23 


12 


1800 


1 


9 


0.44 


13 


00 00 


2 


13 


0.39 


13 


03 00 


9 


13 


0.80 


13 


06 00 


0 


13 


0.80 


13 


09 00 


1 


13 


0.24 


13 


12 00 


1 


13 


0.30 


13 


15 00 


0 


13 


0.31 


13 


18 00 


2 


13 


0.30 


13 


2100 


3 


13 


0.51 


14 


0000 


2 


15 


0.07 


14 


03 00 


2 


15 


0.13 


14 


06 00 


2 


15 


0.24 


14 


09 00 


2 


15 


0.26 


14 


21 00 


3 


15 


0.26 


15 


00 00 


0 


4 


0. 1 6 


15 


03 00 


0 


4 


0.07 


15 


06 00 


0 


4 


0.0 5 


15 


09 00 


1 


4 


0.09 


15 


12 00 


0 


4 


0.24 


15 


15 00 


0 


4 


0. 15 


15 


1800 


1 


4 


0. 19 


15 


21 00 


1 


4 


0.30 


16 


00 00 


1 


8 


0. 1 1 


16 


03 00 


3 


8 


0.48 


16 


/ 06 00 


3 


3 


0.09 


16 


0900 


2 


8 


0.12 


16 


12 00 


0 


8 


0.19 


16 


15 00 


0 


3 


0.24 


16 


1800 


1 


8 


0.1 2 


16 


21 00 


2 


8 


0.30 


17 


00 00 


1 


5 


0.11 


17 


03 00 


3 


* 


0. 14 


17 


06 00 


1 


5 


0.10 


17 


09 00 


0 


5 


0.04 


17 


12 00 


0 


5 


0.08 


17 


1500 


1 


5 


0.10 


17 


1800 


1 


5 


0.07 


13 


00 00 


0 


2 


0.07 


13 


09 00 


0 


2 


0. 1 8 


13 


12 00 


0 


2 


0. 10 


13 


1500 


0 


2 


0.10 


13 


21 00 


0 


2 


0.07 



(after Ochadlick) 
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TABLE VI 



Observed AOTEC Data and Indices, April 19-24, 1976. 



Date 


Time(Z) K- 


San 


Juaa Ap 


Amplitud' 


19 


00 00 


1 


5 


' 0.26 


19 


0300 


1 


5 


0.16 


19 


06 00 


1 


5 


0.05 


19 


09 00 


1 


5 


0.07 


19 


1200 


0 


5 


0.05 


19 


1500 


0 


5 


0.10 


19 


18 00 


1 


5 


0. 1 1 


20 


0300 


1 


2 


0.05 


20 


06 00 


0 


2 


0.04 


20 


0900 


0 


2 


0.08 


20 


1200 


0 


2 


0.03 


20 


15 00 


0 


2 


0.06 


20 


1800 


1 


2 


0.07 


20 


21 00 


0 


2 


0.07 


21 


0000 


0 


5 


0.12 


21 


03 00 


0 


5 


0.04 


21 


0600 


1 


5 


0. 1 1 


21 


09 00 


0 


5 


0.14 


21 


12 00 


0 


5 


0.06 


21 


18 00 


1 


5 


0.11 


22 


00 00 


3 


19 


0.10 


22 


06 00 


4 


19 


0.32 


22 


09 00 


3 


19 


0.21 


22 


12 00 


1 


1 9 


0.36 


22 


21 00 


2 


19 


0.45 


23 


00 00 


2 


8 


0.12 


23 


03 00 


? 


8 


0.17 


23 


06 00 


1 


3 


0.05 


23 


0900 


0 


8 


0.15 


23 


1200 


1 


8 


0.08 


23 


1500 


1 


3 


0.17 


23 


1800 


1 


8 


0.09 


23 


21 00 


2 


8 


0.17 


24 


03 00 


2 


16 


0.33 


24 


06 00 


3 


16 


0.21 


24 


0900 


2 


16 


0.03 


24 


1200 


1 


16 


0. 14 








(after Ochadlic 


The 


largest peak 


-to- 


■peak fluct 


uations o 


ASQ-81 


magnetomete r 


in 


a three 


hour peri' 



compared to the K-San Juan index for that period. Data 

points with fluctuations greater than 1.0 nT were deleted as 
the observatory did not record that information. The sample 
coefficient of linear correlation was 0.34 with the .95 
confidence interval for the actual correlation coefficient 
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being from 0.14 to to 0.52. Sample size was 86. This indi- 
cates that there is at most a weak correlation between the 
observed data and the K-San Juan index. A much greater 
correlation would be required for the K index to be of any 
significant value for HAD operational use. 

The coefficient of correlation has values from 
-1.0 to 1.0. A value of -1 . 0 or 1.0 indicates perfect 
negative or positive correlation, respectively. A value of 
zero signifies no correlation at all. 

A similar correlation analysis against the Ap 
index was conducted. The sample correlation coefficient was 
0.51 with the .95 confidence interval for The actual coeffi- 
cient of correlation being from 0.32 to 0.65. Again, this 
signifies than only a weak correlation was observed. 

The weak correlation between observed geomag- 
netic activity and the K-San Juan index suggests that This 
index would not be very useful in describina HAD geomagnetic 
noise, as this K index reflected activity similar to that in 
the band of interest only about one-third of the time. The 
stonger correlation of the Ap index indicates that it 
reflected activity similar to HAD geomagnetic noise about 
one-half of the time. This is still not a very good indica- 
tion of what is going on in the HAD band pass. 

c. Power Spectral Density Evaluation 

As part of ongoing research at the Naval 
Postgraduate School, geomagnetic activity data in the range 
of the HAD band pass has been collected and analyzed in the 
form of power spectral density (PSD) curves. 

A measure of HAD band activity has been devel- 
oped by integrating under the PSD curve and using compromise 
band pass limits of 0.05 and 1.0 Hz. This type of index is 
discussed in depth later. 



Two sets of data were analyzed. One set was 
taken using an induction coil to measure the fluctuations in 
one direction during day, 1980 [Ref. 39] and August to 
October, 1981 [Ref. 40]. The other data sat was collected 
using a Cesium vapor toatal field magnetometer during the 
July to October, 1980 period [Ref. 41]. 

Correlation analysis of the single-coil system 
data yielded a sample correlation coefficient of 0.401 for 
the K -Fredericksburg index, 0.180 for the A-Fredericksburg 
index and 0.155 for the Ap index. Sample size was 9. The 
.95 confidence intervals for the correlation coefficient for 
K-?red, A-Fred, and Ap ware -.36 to .85, -.55 to .76, and 

-.56 to .75, respectively. The data for this test is 
presented in Table VI I. 

TABLE VII 

Single Coil RdS Noise Data and Indices 











RdS Noise 






Dat 


6 


Time { Z) 


Amol if ude (nT) K-Fred 


A-Fred Ap 


1 


day 


80 


1800 


'0.10 3 


10 8 


2 1 


Auer 


81 


2045 


0.05 2 


12 14 


25 


Aug 


81 


0235 


0. 12 3 


15 16 


25 


Auq 


31 


0745 


0. 19 3 


15 16 


25 


Audr 


8 1 


1318 


0.05 3 


15 16 


4 


Oct 


8 1 


1830 


0.06 3 


17 17 


5 


Oct 


81 


0155 


0. 12 1 


7 7 


5 


Cct 


81 


0630 


0.09 3 


7 7 


5 


Oct 


8 1 


1530 


0.02 0 


7 7 








The Cesium 


vapor magnetometer 


data yielded 


sample ccef 


ficients of 


correlation of 0.552 


, 0.374, and 


0.4 44 


fc: 


r the K-Fred, A 


-Fred, and Ap indices. 


respectively. 


The s 


ample 


size was 14. 


The .95 confidence 


intervals were 


for K 


-Fred 


from .03 to 


.84, A-Fred from -.20 


to .76, and for 


Ap fr 


om 


“ • 


12 to .79. 


The Cs vapor magnet 


ouster data is 



presented in Table VIII. 
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TABL2 VIII 



Cs Vapor BflS Noise Data and Indices, Jul-Oct , 1980 







RMS Noise 






Date 


Time (Z) Ampl itude (nT) 


K-Fred 


A- Fred Ap 


2 5 Jul 


1800-1830 


0. 10 


6 


31 46 


3 0 Jul 


1450-1520 


0.06 


2 


7 8 


30 Jul 


1520-1551 


0.07 


2 


7 8 


3 1 Jul 


14 50-1520 


0.05 


1 


6 6 


31 Jul 


1520-1555 


0.04 


1 


6 6 


16 Aug 


0030-0100 


0. 05 


2 


19 26 


16 A uq 


0900-0930 


0.04 


2 


19 26 


18 Oct 


02 30-0300 


0.06 


3 


17 21 


1 8 Oct 


0300-0330 


0.06 


4 


17 .21 


18 Cct 


0630-0700 


0. 04 


4 


17 21 


18 Oct 


1030-1 ICO 


0. 04 


3 


17 21 


18 Oct 


1 100-1 130 


0. 05 


3 


17 21 


18 Oct 


1 2 30- 1 30 C 


0.04 


3 


17 21 


18 Oct 


1300-1330 


0.04 


3 


17 21 




Although the 


sample sizes 


used were too small to 


draw any meaningful conclusions, there is little evidence to 


suggest tha 


t any of the 


K-Fred, A-?r 


ed , or 


Ap indices is a 


very accura 


he measure of 


ge omagnetic 


noise in the ASQ-91 MAD 


band pass. 










a. 


Correlation 


Con clus ions 








Although some weak co 


rr elation does exist 


between the 


S indices. 


the A-Fred i 


naex, i 


:hs Ap index and 


geo magnetic 


noise in the 


MAD band pas 


s, this correlation is 


incidental 


and indirect. 


being the 


result of a correlation 


between th 


e activity i 


n the MAD band and in the lower 



frequency activity that influences the K and A indices. 
These indices are not directly influenced by activity in the 
MAD band. The correlation that does exist does not appear to 
be sufficiently high to enable these indices to yield accu- 
rate indications of the actual MAD band activity. The use 
of these indices for anything except the roughest qualita- 
tive estimation of active in the MAD band pass is not 
recommended. 
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D. PROPOSED GEOMAGNETIC INDICES FOR MAD 

Overall geomagnetic acitivity is analyzed in both the 
time and frequency domains. Geomagnetic noise indices for 
the MAD band pass could be developed in either of these 
domains. 

The spatial coherence of MAD geomagnetic noise has not 
yet been adequately determined. This information would be 
necessary in order to determine the number and location of 
mini- observatories for an operational MAD noise index. 

1 . Time Seri es A nalysi s 

One way to develop an index of MAD geomagnetic 
activity would be to establish mini-observatories near bases 
from which MAD operations are conducted. These observato- 
ries would use ASQ-81 magnetometers or different magnetome- 
ters with ASQ-81 filter networks, and could in real time 
record the geomagnetic noise in the MAD band. A measure 
such as the maximum peak-to-peak (or possibly the average 
peak-to-peak) noise in a given time period could then be 
disseminated to flight craws operating in the area covered 
by that index. Obviously, the spatial coherence of MAD band 
activity is important in making such a system work. This 
type of mini-observatory has been suggested by References 31 
and 3 5. 

2. Frequen cy Dom ain Index 

Present fleet procedures examine MAD noise such as 
system and manuever noise in terms of the amplitude of the 
fluctuation [Ref. 30]. An index of geomagnetic noise in th° 
MAD band pass would therefore be of greatest usefulness to 
the fleet operator if it were in units of the amplitude of 
the signal as seen by the MAD equipment. 
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The method proposed to derive a MAD noise index in 
the frequency domain begins with obtaining the power spec- 
tral density of the activity in the MAD band by Fourier 
analysis of the time series data input from the 
magnetometer. 

By intergrating under the PSD curve over the limits 
of the MAD bandpass, a value in units of amplitude 2 will 
result. Taking the square root of ^his value will yield an 
RMS amplitude. Equation 4-1 represents the derivation of 
this index. 



MAD Index 




PSD (f) 




1/2 



(4-1) 



where MAD 
f 
u 
1 

and PS D ( f ) 



Index is in nT (gammas) 

= frequency (Hz) 

= upper bandpass limit 
= lower bandpass limit 
= Dower soectral density 
(r.T 2 /hz) 



The characteristics of the MAD filter (the filter 
not being an ideal bandpass filter) could be applied prior 
to the integration. The integration itself could be done by 
either a point-bv- point numerical integration or by first 
modelling the PSD curve by polynomial curve fitting and then 
integrating the polynomial over the range of the band pass. 
It is anticipated that this could be done in close no real 
time by a digital computer, possibly by a desk top computer 
such as the HP9845. 

The type of sensor utilized could be the ASQ-81 
magnetometer, other total field magnetometer, or possibly an 
orthogonal 3-coil system whose signals can be combined to 
yield the projection on the total field vector of the fluc- 
tuations. A single coil syszem oriented in the direction of 
the earth's magnetic field vector could also be used. 
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a three coil system which is used to yield an RMS 
amplitude is currently in research use at the Naval 
Postgraduate School. 

3. P redict io ns of Geomagnetic Ac try it v 

The proposed indices discussed above are intended to 
be real time measures of the geomagnetic noise in the MAD 
bandpass. Whether or not such activity can be predicted 
ahead of time needs to be looked info. 

While there is no model for the background component 
of geomagnetic noise, work has been done on estimating the 
future activity of micropui sations, notably by Praser-Saith 
in the case of Pci pulsations [Ref. 42, 43]. 3y extending 

the prediction technique for Pci pulsations to the Pc2, Pc3, 
and Pil pulsations, the occurrence of geomagnetic micropul- 
sation activity in the MAD band might be predicted. 
Combining this prediction with real-time solar flare infor- 
mation should give the capability to disseminate real-time 
and estimated future MAD index values to fleet users. 
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conclusion that 



Previous studies have led to the conclusion that the 
currently used geomagnetic indices ace not accurate measures 
of geomagnetic noise in the MAD band pass. Experimental 
equipment has been utilized, and computer software written 
in order to confirm this conclusion, and to develop a 
replacement means of evaluating MAD geomagnetic noise. 



A. EQUIPMENT CONFIGURATION 

Experimental equipment, acquired as part of the Naval 
Postgraduate School gecmagnetics research program, has been 
utilized in the effort to develop a usable MAD index. This 
equipment is in use in other projects of the gecmagnetics 
research group. The sensors and associated equipment are set 
up for remote site operation with system monitoring and data 
analysis located at the Naval Postrgraduat e School. 
Descriptions of the data collection system and data analysis 
system fellow below. 



1 . Ba^a Collection Sy stem 

The data acquisition system illustrated in Figure 
5.1 reveals the following major components: 

-ceil antenna sensors (3) 

-preamplifiers (3) 

-signal conditioners (amplifiers) (3) 

-pulse code modulation system (1) 

-radio transmitter (1) 

-radio receiver (1) 

-instrumentation tape recorder (1) 
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PREAMP — — i AMP 



r 




Figure 5.1 Data Collection System 
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PWR 




a. Coil Antenna Sar.sors 

Three coil sensors are used in this system. 
Each sensor is a self-supporting, continuously wound, non- 
center-tapped coil antenna manufactured by Sima Engineering, 
Palo Alto, California, from about 5460 turns of 18 gauge 
copper magnet wire. The coils weigh approximately 50 kg each 
with dimensions as depicted in Figure 5.2. The dimensions 
of the sensor are constrained by the dimensions of the 
largest glass sphere that is commercially available. These 
spheres are used to enclose the coils during underwater 
experiments. The coil resistance is 120 ohms and its self- 
inductance is approximately 9.31 henries. The three coils 
are mounted orthogonaly on a nonmagnetic frame (Figure 5.3). 




Figure 5.2 Sensor Dimensions 
b. Preamplifier 

The preamplifier used was the model 1 3 — 1 0 A lew 
noise ELF amplifier manufactured by Dr. Alan Phillips of SRI 
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Figure 5.3 Sensor Mounting Block 



International. The final stage of the amplifier contains an 
active low-pass filter which provides a sharp curoff for 
frequencies above 20 Hz. The overall preamplifier gain for 
inputs of less rhan 2.5 millivolts is 60 dB. 



c. Signal Conditioners 

The signal conditioners receive the analog 

signals from the coil preamplifiers, amplify them on the 

order of 30 dB, and limits signals with peak amplitudes of 
7.5 volts from entering the pulse code modulation system. 

d. Pulse Code Modulation (PCM) System 

The pulse cods modulation system chosen for use 
is one designed and manufactured by Dr. Robert Lowe, 
Lowecom, Inc. The system features 15 channel analog input 
capability and offers selectable sampling rages of 2, 4, 8, 

16, 32, 64, and 128 samples per second. By appropriately 

jumpering the analog input pins, the sampling rate may be 
increased by a factor of 5. 

The PCM system incorporates a crystal oscil- 
lator, and associated CMOS integrated circuitry tc develop 
the clocking pulses, and a 16 channel CMOS analog multi- 
plexer, a 16 channel, 12 bit CMOS analog to digital 
converter and associated circuitry to provide the pulse 
coding. The crystal clock oscillator operating at a 
frequency of 24.576 kHz produces a square wave output with a 
loss rate of 1 bit in 10 6 . The clock pulses gate the analog 
multiplexer, analog to digital converter and associated 
follow cn circuitry that form the pulse code words. The 
basic output is a 3i-phased pulse coded signal. 

The data is organized in frames. Each frame is 
headed by a sync code word which is followed sequentially by 
the pulse coded samples from PCM channels 1 through 15. The 
sync code word is a pulse coded digital word with a decimal 
value between 0 and 4095. This word is preselected and 
hardwired on the circuit board. This code word is essential 
to the decoding process. 
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In the initial operation of this system, a 
sampling rate of 32 samples per second was utilized. Only 
one sample per coil per frame was analyzed, making use of 
PCM channels 2, 3, and 4 only. 

e. Transmission and Recording 

After the data has been PCM encoded, it is 
transmitted by a VHF radio link back to a receiver located 
at the Naval Postgraduate School, where it is currently 
recorded on an instrumentation tape recorder for iater 
analysis. 

2. Data Anal ysis Equi p me nt 

Currently, the recorded PCM data is played back into 
a PCM decoder and associated equipment which generates a 
nine-track 800 BPI computer tape containing the decoded 
sensor data. This computer tape is then input into the 

Naval Postgraduate School IBM 3033 computer for analysis. It 
is in the computer software that the senscr system transfer 
function is applied, spectral analysis performed, and the 
MAD index generated. 

3. DATA ANALYSIS SOFTWARE 

As was noted earlier, a mainframe computer was utilized 
to perform the spectral analysis, apply the transfer func- 
tion, convert to power spectral density, plot the PSD and 
generate the RMS noise amplitude MAD index (using equation 
4-1). This program is written in FORTRAN IV and is 
discussed in brief below. A copy of the program can be 
found in Appendix A. 

The main program is divided into sections which perform 
the following functions: 

-Data input 

-Fourier analysis of time series data 
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-Application of system transfer function 

-Projection of field components onto total field vector 

-Data averaging 

-Curve fitting and calculation of MAD index 
-Calculation and plotting of power spectral density 

1 . Data In p ut 

Data input is accomplished with the aid of a subrou- 
tine package supplied by Dr. Tim Stanton of the Naval 
Postgraduate school Department of Oceanography. His subrou- 
tine (called 'SUBROUTINE SD') serves as a FORTRAN 'READ' 
statement, taking the PCM data off the computer tape and 
converting it info integer format with a value between 0 and 
4095. The data input section of the main program takes this 
integer value and converts it to a 'REAL' number and normal- 
izes it to represent a voltage value between -5.0 and +5.0 
volts. This section also sorts the input data matching the 
PCM channel to the data array representing the appropriate 
coil. 



2 • Fouri er A nalysis 

The time series data is next converted to the 
frequency domain by utilizing a subroutine (called 'FOURT') 
which performs a Fast Fourier Transform (FFT) . The subrou- 
tine is one available to users at the Naval Postgraduate 
School and utilizes che Cooley-Tukey FFT algorithm. Further 
information about this subrourin® can be found in the 
program listing in Appendix A. 

3 . Ap plicati on of Tra nsfe r ? unction and Total Field 

Project ion 

The next section of cede applies the system transfer 
functions for the three coils to the frequency domain aa^a. 
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Ths transfer functions are given as straight line segments 
which were found by least- squares approximation. The data 
enters this section in amplitude units of volts, and is 
converted into nanoTeslas by the transfer function. 

Following the application of the system transfer 
function to the coil data, the program next calculates the 
projection of this data onto the earth's magnetic field 
vector (total field projection). This is done by first 

applying the local magnetic variation (declination) to the 
North-South (X) and East-West (Y) coil information to deter- 
mine the horizontal field component. The local magnetic dip 
(inclination) angle is then used to project the vertical 
(Z-coil) and horizontal field fluctuations onto the total 
field fluctuation. 

4 . Data Aver ag in g 

The previous program, sections exist inside of a 
do- loop which enables the analysis of a long period of data 
without a prohibitive need for storage space. This loop 
includes accumulator arrays for each field component and the 
total field projection. The fluctuation data is converted 
into power prior to storage. This is done by faking the 
magnitude of the fluctuation component, dividing by the 
number of sample points, and then squaring the value. After 
the program passes through the averaging loop for the last 
time, the arithmetic average is taken for each frequency 
point cn the arrays. At his stage the power spectrum is 
multiplied by the sample period to determine power spectral 
density. 

5. M2 Index Calculation 

The next section computes the RMS NAD noise index 
previously discussed. A polynomial curve fit is performed on 
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the total field PSD using an available subroutine ( ' CH3FT ' ) . 
The resulting polynomial is then integrated over the limits 
of the various ASQ-81 band pass settings. 

6 . Plo tting of Power Spectral De ns ity 

Flots of the power spectral density of each of the 
field components and the total field projection are gener- 
ated in the last section of the program. This is dene by 
converting the fluctuation power spectral density to deci- 
bels (dB) referenced to 1 nanolesla 2 per Hertz. A Vsrsatec 
plotting subroutine ('PLOTP') is then called to actually 
generate the plots. 

C. INITIAL SYSTEM OPERATION 

The NPS MAD index system was initially placed into oper- 
ation with the coil sensors located in the La M«=sa Village 
housing area near the Naval Postgraduate School, Monterey, 
California. System checkout was accomplished in June, 1982. 
The full system was placed into operation on 25 July 1982 
and 18 August 1982 in conjunction with similar measurements 
taken on the floor of Monterey Bay. 

The MAD index output and power spectral density plors of 
the total field fluctuation for 25 July 1992, 1237-1406 

local (2037-2206Z) , and 18 August 1982, 0121-0250 local 

(0921-1050Z) and 0507-0636 local (1307-1436Z) are shown in 
Tables IX, X, XI, and Figures 5.4, 5.5, and 5.5 

respectively. 
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TABLE IX 



a AD Noise 2037- 


-22 06 Z , 25 JOL 32, Monterey, CA 


1 Gamma = 


= 1 NanoTesla 


MAD IND2X= 
BANDPASS: 


0. 1083 GAMMAS 
0.04 TO 0.20 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1148 GAMMAS 
0.04 TO 0.40 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1169 GAMMAS 
0.04 TO 0.60 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 0888 GAMMAS 
0.06 TO 0.20 HERTZ 


MAD INDEX* 
BANDPASS: 


0.0966 GAMMAS 
0.05 TO 0.40 HERTZ 


MAD INDEX* 
BAN DPASS: 


0. 0991 GAMMAS 
0.06 TO 0.60 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 0729 GAMMAS 
0.08 TO 0.20 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 0822 GAMMAS 
0.08 TO 0.40 HERTZ 


MAD INDEX* 
3 AN DPASS: 


0. 0352 GAMMAS 
0.08 TO 0.60 HERTZ 


MAD INDEX* 
BANDPASS : 


0. 0598 GAMMAS 
0.10 TO 0.20 HERTZ 


MAD INDEX* 
BANDPASS : 


0. 0708 GAMMAS 
0. 10 TO 0.40 HERTZ 


MAD INDEX* 
3 ANDPASS : 


0.0743 GAMMAS 
0.10 TO 0.60 HERTZ 


MAD INDEX* 
3 AN DP AS S : 


0. 1234 GAMMAS 
0.04 TO 2.00 HERTZ 


MAD INDEX* 
3 ANDPASS : 


0. 1122 GAMMAS 
0.06 TO 2.00 HERTZ 


MAD INDEX* 
BANDPASS : 


0. 1018 GAMMAS 
0.08 TO 2.00 HERTZ 


MAD INDEX* 
3 AN DPASS : 


0. 0923 GAMMAS 
0.10 TO 2.00 HERTZ 
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Figure 5.4 PSD 25 JOL 32, 2037-2206Z, La Mesa Village 



79 



-uo 



TABLE X 



MAD Noise 0921- 


-1050Z , 18 AOG 32, Monterey, CA 


1 Gamma : 


= 1 NanoTesla 


MAD INDEX= 
3 ANDPASS : 


0. 3644 GAMMAS 
0.04 TO 0.20 H BETZ 


MAD INDEX= 
BANDPASS: 


0. 3960 GAMMAS 
0.04 TO 0.40 HERTZ 


MAD INDEX= 
BANDPASS: 


0. 4004 GAMMAS 
0.04 TO 0.60 HERTZ 


MAD INDEX= 
BANDPASS : 


0. 3027 GAMMAS 
0.06 TO 0.20 HERTZ 


MAD I NDEX= 
BANDPASS : 


0. 3401 GAMMAS 
0.06 TO 0.40 HERTZ 


MAD I NDEX= 
3 ANDPASS : 


0. 3453 GAMMAS 
0.06 TO 0.60 HERTZ 


MAD INDEX= 
BANDPASS: 


0.2528 GAMMAS 
0.08 TO 0.20 HERTZ 


MAD I NDEX= 
3ANDPASS : 


0. 2966 GAMMAS 
0.08 TO 0.40 HERTZ 


MAD I NDEX= 
3 AN DP ASS : 


0. 3025 GAMMAS 
0.08 TO 0.60 HERTZ 


MAD INDEX= 
BANDPASS : 


0.2118 GAMMAS 
0. 10 TO 0.20 HERTZ 


MAD INDEX= 
BANDPASS: 


0. 2624 GAMMAS 
0. 10 TO 0.40 HERTZ 


MAD INDSX= 
3 AN DP ASS : 


0. 2691 GAMMAS 
0. 10 TO 0.60 HERTZ 


MAD INDEX= 
3ANDPASS: 


C. 4047 GAMMAS 
0.04 TO 2.00 HERTZ 


MAD INDSX= 
3 ANDPASS : 


0. 3679 GAMMAS 
0.05 TO 2.00 HERTZ 


MAD INDEX= 
BANDPASS: 


0.3339 GAMMAS 
0.08 TO 2.00 HERTZ 


MAD I NDEX= 
BANDPASS: 


0. 3026 GAMMAS 
0.10 TO 2.00 HERTZ 
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Figure 5.5 PSD 18 AOG 82, 0921-1050Z, La Mesa Village 



81 



-60 



TABLE XI 



HAD Noise 1307- 


-14 36 Z , 18 AUG 32, Monterey, CA 


1 Gamma = 


= 1 NanoTesla 


MAD INDSX= 
BANDPASS: 


0. 3188 GAMMAS 
0.04 TO 0.20 HERTZ 


MAD I NDEX* 
BANDPASS: 


0.3246 GAMMAS 
0.04 TO 0.40 HERTZ 


MAD INDEX* 
BANDPASS : 


0. 3257 GAMMAS 
0.04 TO 0.60 HERTZ 


MAD TNDEX= 
BANDPASS: 


0. 2297 GAMMAS 
0.06 TO 0.20 HERTZ 


MAD INDEX= 
BANDPASS: 


0. 2376 GAMMAS 
0.06 TO 0.40 HERTZ 


MAD INDSX= 
BANDPASS : 


0. 2405 GAMMAS 
0.06 TO 0.60 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1585 GAMMAS 
0.08 TO 0.20 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1693 GAMMAS 
0.08 TO 0.40 HERTZ 


MAD INDEX* 
BAN DPASS: 


0. 1738 GAMMAS 
0.08 TO 0.60 HERTZ 


MAD INDEX* 
3 AN DPASS : 


0. 1031 GAMMAS 
0.10 TC 0.20 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1199 GAMMAS 
0.10 TO 0.40 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 1254 GAMMAS 
0.10 TO 0.60 HERTZ 


MAD INDEX* 
3 AN DP ASS : 


0. 3504 GAMMAS 
0.04 TO 2.00 HERTZ 


MAD INDEX* 
BAN DPASS: 


0.3137 GAMMAS 
0.06 TO 2.00 HERTZ 


MAD INDEX* 
BANDPASS: 


0. 2801 GAMMAS 
0.03 TO 2.00 HERTZ 


MAD INDEX* 
BANDPASS : 


0. 2495 GAMMAS 
0. 10 TO 2.00 H3P.TZ 
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VI 



• CONCLUSION AND RECOMME NDATI ONS 

A. CONCLUSION 

Analysis of available information regarding current 
geomagnetic indices and the actual level cf geomagnetic 
noise in the MAD bandpass indicates than the currently used 
indices, the K and Ap indices, are not valid for MAD opera- 
tions. It is therefore desirable to derive a new index 
which more accurately represents the geomagnetic noise at 
frequencies of interest in MAD operations. 

3. RECOMMENDATIONS 

A tentative index was developed, tested, and sample data 
obtained. The data analysis for the preliminary system was 
accomplished using a mainframe computer. It is recommended 
that further work include improving the index and the 
setting up of a system to on-line decode the incoming data, 
and utilizing a desk top computer such as the HP9845 to 
enable real time determination of the MAD noise index. 

The spatial coherence of MAD geomagnetic noise should 
also be investigated, with possible application to noise 
cancellation. Additionally, the feasibility of the prior 
estimation of geomagnetic noise should be evaluated, 
possibly by using the prediction methods proposed by 
Eraser-Smith [Ref. 42, 43]. 
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MAD INDEX DATA ANALYSIS SOFTWARE 
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